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Abstract. The electricity market is a very peculiar market due to the large variety 
of phenomena that can affect the spot price. However, this market still shows 
many typical features of other speculative (commodity) markets like, for instance, 
data clustering and mean reversion. We apply the diffusion entropy analysis (DEA) 
to the Nordic spot electricity market (Nord Pool). We study the waiting time 
statistics between consecutive spot price spikes and find it to show anomalous scaling 
characterized by a decaying power-law. The exponent observed in data follows a quite 
robust relationship with the one implied by the DEA analysis. We also in terms of the 
DEA revisit topics like clustering, mean-reversion and periodicities. We finally propose 
a GARCH inspired model but for the price itself. Models in the context of stochastic 
volatility processes appear under this scope to have a feasible description. 
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1. Introduction 

During tlie last years, there has been an increasing number of contributions from the 
physics community to the study of economic systems. Energy spot prices, that result 
from the deregulation of the power sector, are no exception. Weron and Przybylowicz 
and Weron P] deal with the Hurst exponent (or R/S) analysis |H:, the detrended 
fluctuation analysis (DFA) |lj and periodogram regression methods. These techniques 
were used to verify and quantify a claim already stated by financial mathematics: 
electricity spot prices are mean reverting OIHl- This means that they suffer a strong 
restoring force driving the price toward a certain normal ("fundamental") level. Using 
the language of physics, one says that prices are anti-persistent, or equivalently, that 
the price increments are negatively correlated. 

Recently, also the Average Wavelet Coefficient (AWC) method [Zj has been apphed 
to spot prices 0. This method shows its potential in particular when dealing with multi- 
scale time series. Due to its separation of scale property, the presence of one scaling 
regime covering a given time range is not hampered by the presence of another one. This 
is not always the case for many other method, like the DFA, where one scaling behavior 
can "spill-over" to the next one and even fully destroy the scaling property of the latter 
one [H]. In terms of power markets, this is highly relevant since the statistical behavior 
of the price on an intra-day scale is mainly determined by the consumption patterns, 
and does not (on this scale) show the characteristic mean reversion character that can 
be observed on the day-to-day scale, and above |H] . The lack of possibility in separating 
the various scaling regimes, with the technique used, was the main motivation for e.g. 
Refs. 12] analyzing mean daily data, which show only one scaling region, instead of 
the original hourly data. 

All the analyzing techniques mentioned above relay on the scaling of some kind 
of fluctuation measure, say, the standard deviation or variance, as function of the 
window size. Such approaches will only measure the correct correlation scaling exponent, 
i.e. the Hurst exponent, if the underlying time series is consistent with Gaussian 
statistics 0. This is for instance the case for the celebrated fractional Brownian 
motion [TOI. For correlated non-Gaussian increments:}: the associated exponent will 
partly receive contributions from the correlations as well as the non-Gaussian character, 
making it difficult (or impossible) to separate the two ^9J. This latter situation is faced 
in e.g. ordinary Levy walks and in its fractional equivalent fO]- Recently, a method 
was proposed that in a reliable way can determine and separate the contribution to 
the scaling from both correlations and non-Gaussian statistics. This method is based 
on the thermodynamics of the time series and known as the diffusion entropy analysis 

(DBA) [nma. 

In this paper we apply the DEA for the study of the electricity spot prices from 
the Nordic Power exchange (Nord Pool). One of our main aims is to investigate the 
statistics of the waiting times between consecutive spot price spikes of deregulated 

I It is here meant that the tails of the distribution are fatter than those of the Gaussian distribution. 



Entropy of the Nordic electricity market 



3 



electricity markets, and to uncover if they show anomalous non-Poissonian statistics. 
The DEA is specially suited for intermittent signals, i.e., for time series where bursts 
of activity are separated from periods of quiescent and regular behavior. The technique 
has been designed to study the time distribution of some markers (or events) defined 
along the time series and thus discover whether these events satisfy the independence 
condition ^TJ ^]. Other objectives are the study of the antipersitency behaviour in 
data via the DEA technique, the observation of existing periodicities in data in many 
different ways and finally the proposal of a GARCH model showing properties similar 
to those of real data. 

This paper is organized as follows. We start in Sec. |21 by briefly discussing the 
spot market and the data set to be analyzed in this work. Then we briefly review the 
DEA technique (Sec. Ej). Section HI studies the statistics of the most important price 
movements and infer some properties on the inter-event peak probability. Section El 
presents the results obtained by the DEA technique on tick-by-tick spot data and 
without filtering the most relevant price changes. Furthermore, we propose in Sec. IHl 
a GARCH model [131 CH 113 UHl HI] for the spot electricity price and try to obtain 
consistent entropies for the tick-by-tick data and the spike filtering. Conclusions are left 
for Sec. [7| 

2. Nord Pool — the Nordic power exchange 

The Nordic commodity market for electricity is known as Nord Pool ^H]- It was 
established in 1992 as a consequence of the Norwegian energy act of 1991 that formally 
paved the way for the deregulation of the electricity sector of that country. At this time it 
was a Norwegian market, but in 1996 and 1998 Sweden and Finland joined, respectively. 
With the dawn of the new millennium (2000), Denmark decided to become member as 
well. 

Nord Pool was the world's first international power exchange. In this market, 
participants from outside the Nordic region are allowed to participate on equal terms 
with "local" participants. To participate in the spot market it is required that the 
participants must have an access to a grid connection enabling them to deliver or take 
out power from the main grid. For this reason, the spot market is often also called the 
physical market. As for today, the physical market has a few hundreds of participants. 
More than one third of the total power consumption in the Nordic region is traded in 
this market, and the fraction has steadily been increasing since the inception of the 
exchange in the early 1990s. In addition to the physical market, there is also a financial 
market. Here, power derivatives, like forwards, futures, and options are being traded. 
This market presently has about four hundred participants. For each of these two 
markets, about ten nationalities are being represented among the market participants. 

Nord Pool is an auction based market place where one trades power contracts for 
physical delivery within the next day. It is known as a spot market. However, in a strict 
sense this notion is not precise since formally it is a day- ahead (24- hours) forward market. 
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Figure 1. The Nord Pool system spot price from May, 1992 till the end of 2004. 
Several (seasonal, weekly, and daily) periodicities can be observed in this data set 
below the spiky randomness. The inset shows the price variations over a (randomly 
chosen) weekly period. 

What is traded are one-hour-long physical power contacts, and the minimum contract 
size is O.lMWh. By noon (12:00 hours) every day, the market participants submit their 
(bid and ask) offers (including prices and volumes) to the market administrator (Nord 
Pool). The offers are submitted for each of the individual 24 hours of the next day 
starting at 1:00 hours. After the submission deadline (for the next day), Nord Pool 
proceeds by preparing (for ever hour) cumulative volume distributions (purchase and 
sale curves) vs. price (p) for both bid {Vb{p)) and ask {Va{p)) offers. Since there in the 
electricity market must be a balance between production and consumption, the so-called 
system (spot) price, P{t), for that particular hour (t) is determined as the price where 
Va(-P) = Vb{P)- This is called the market cross, or equilibrium point. Trading based on 
this method is called equilibrium trading, auction trading or simultaneous price setting. 
If the data do not define an equilibrium point, no transactions will take place for that 
hour, and no system spot price will therefore be determined. So far, to our knowledge, 
this has never happened at Nord Pool. 

After having determined the system price, P{t), for a given next-day hour, Nord 
Pool looks for potential bottlenecks (grid congestions) in the power transmission grid 
that might result from this system price. If no bottlenecks are found, the system price 
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will represent the spot price for the whole Nord Pool area for that given hour. On the 
other hand, if potential grid congestions may result from the bidding, so-called area 
(spot) prices, that are different from the system price, will have to be created. The idea 
behind the introduction of area prices is to adjust electricity prices within a geographical 
area in order to favor local trading to such a degree that the limited capacity of the 
transmission grid is not exceeded. How the area prices are being determined within 
Nord Pool differs between, say, Sweden and Norway, and we will not discuss it further 
here (see e.g. Ref. [IH] for details). 

In this work, we will analyze the hourly Nord Pool system spot prices for the 
period from (Monday) May 4th, 1992 till the end of Friday December 31st, 2004; in 
total 110,987 data points (Fig. H]). In Fig. |21we present the corresponding normalized 
hourly returns (to be defined in Eq. ^ below). The reader should notice that high levels 
of return are possible, and typical, for electricity markets. From Fig. Q one should also 
observe that the price process shows several periodicities. Those are mainly attributed 
to consumption patterns and are daily, weekly and seasonal in character. They have 
been reported and studied previously by several research groups |H1 UHl UH 1201 HI] 
Superimposed onto this deterministic periodic trend, is a random component showing 
strong variability with pronounced spikes and data clustering. The DEA technique 
allows us to deal with this randomness by somewhat ignoring the periodic signal. 

3. The Diffusion Entropy Analysis 

The DEA technique is a statistical method for measuring scaling exponents of time 
series by utilizing their thermodynamical properties ^Tj ^]. This is achieved by (i) 
converting the time series into some kind of probability density function (pdf), p{x,t), 
where the variable x is related to the fluctuating time series and t denotes the time (or 
time interval), and (ii) therefrom calculating the related Shannon (information) entropy 
S{t) (to be defined below), from which the scaling properties (if any) of the pdf can be 
deduced (c/. Eqs. ((H) and below). 

Various pdf's of the form p{x,t), obtained from the underlying time series, can be 
used together with the DEA technique. For instance, p{x, t) can be the probability of 
being at position x at (diffusion) time t given that one was at x = at time t = 0, or x 
may be a "marker" related to the number of events (like values of the time-series being 
above a given threshold etc.) occurring within a time interval of length t. Having said 
that, it should be mentioned that generally x should be a zero-mean variable in order to 
avoid artifacts in the application of the DEA fH^]. We will in the preceding section 
describe in more details how one may construct a suitable pdf to be used in the DEA. 

Under the assumption that the time series is stationary and scale invariant, one 
has that the pdf can be written as 



where S denotes the so-called scaling exponent (that one wants to determine), and F is 




(1) 
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Figure 2. The normalized Nord Pool hourly {At = 1 hour) (logarithmic) returns, 
fAt{t)/<^At (c/- Eq. lO), of the data depicted in Fig. ^ The normalization is done 
according to the average hourly volatility (TAt- Notice that normalized return values 
higher than 10 is rather frequent and typical for this (and many other) electricity 
markets. 



some positive and integrable function depending on the specificities of the pdf. Then 
the Shannon entropy, that in general is defined as 

/oo 
dxp(x,t) ln[p(x,t)], (2) 
-oo 

wiU take on the form 

S{t)=A + 5\nt, (3) 

where A is given by an expression similar to Eq. ((2)), but with F{y) substituted for 
p{x,t). This last transition is easily demonstrated by substituting Eq. (jT)) into Eq. (j2I) 
and making a change of variable to y = x/t^. 

The important lesson that should be taken from Eq. ^ is that the scaling exponent, 
6, can be obtained as the slope of the entropy vs. \nt curve. This is a reliable method 
13 HH nil for measuring the scaling exponent of the probability density (H} . 

4. Spikes and data clustering 

One of the most pronounced features of the electricity spot price process of Fig. [T] is its 
spiky nature. Within one, or only a few hours, the spot price can increased manifold. 
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Figure 3. The frequency distributions of the 100 largest hourly return events vs. 
month of year and hour of the day (inset) based on the data depicted in Fig. El This 
result clearly indicates that the dramatic price changes are not evenly distributed in 
time, but instead occur during periods of high consumption. 



One of the most extreme example of such a situation occurred on February 5, 2001, when 
the spot price reached an all-time-high of 1951.76 NOK/MWh. Three hours earlier, the 
spot price had been at about 200 NOK/MWh, a more "normal" level for that time of 
year. 

The dramatic price variations take place during periods of high consumption, 
which for the Nordic area means during the winter season as illustrated by Fig. IHl 
One way to explain the appearance of these dramatic price changes is by the so- 
called stacking or inventory models j22]: Extra demand for electric power is normally 
filled by the cheapest available energy source. During low consumption periods, the 
daily consumption profiles, say, will not influence the price level in any dramatic way, 
since there is plenty of base energy generation capacity available. However, when the 
consumption already is high, only a minor increase in the demand can have dramatic 
consequences for the electricity prices. The higher the demand, the more costly will 
it typically be to produce an extra unit of energy, and the more sensitivity will the 
price of electric power be to minor changes in the consumption. If one during such 
high-consumption periods experiences loss of less costly generating capacity, e.g. due to 
technical problems, extreme situations like the one reported above can be experienced. 
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This sensitivity of the spot price to the level of consumption implies seasonal volatility 
clustering with typical low and high volatility periods during the summer and winter 
seasons, respectively [THj . 

The fact that the extreme price spikes are more likely to occur during the winter 
season than during the remaining parts the year, has the immediate consequence that 
time intervals between (consecutive) spikes are not expected to constitute a Poisson 
process. Hence, the time intervals are expected to scale anomalously according to 
Feller [23], and hence the scaling exponent will be different from 1/2. This waiting 
time statistic (between spikes) has previously been studied in e.g. Ref. P for solar 
flares or earthquakes. It is worth mentioning that although the presence of seasonalities 
in spike statistics is quite evident looking at Fig. El this does not mean that the spike 
events follow a somewhat deterministic rule. As we will show later, time between spikes 
T appears to be described by a power law probabihty distribution. The periodicities 
affect the statistics of very particular values of r but not the whole probability density 
function. 

It has in the past been demonstrated ^2] that if the "spike" events are independent 
and the time interval between them distributed (at least asymptotically) according to 
a power law of the form {(3 > 0) 

0(r) ~ I/t", (4) 

then an analytic relationship between the (3 and the scaling exponent 6 do exist. For 
instance, if one constructs an auxiliary process by introducing steps, say, of magnitude 
1, whenever the originally set has absolute returns larger than a given threshold (called 
a marker below), the relationship reads J2] 

f /3 - 1, if 1< ^ < 2; 

5spikes= I l/iP-l), if 2</?<3; (5) 
[l/2, if P>3. 

Under this scope, the present section investigates the predictive power of the DEA 
analysis applied to the spike statistics of the Nordic electricity market. From the hourly 
system spot prices data, P{t), one defines the logarithmic returns§ over the time horizon. 
At, as 

and in the following it will be assumed that At = 1 hour. We intend to investigate 
the properties and time distribution of those returns being larger than a certain lower 
threshold tq. To this end, we define a time-dependent marker position variable by 

«t) = |'- 'f > '■»• (7) 

I 0, otherwise. 

§ The reader should notice that due to the large price variations (relative to the price level itself) 
that can occur for electricity spot prices, relative and logarithmic returns are not necessarily even close 
to being similar. This is in sharp contrast to, say, stock markets where the two types of returns are 
approximately equal in most cases. 
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Table 1. Regressions for the diffusion entropy results given in Fig. ^ assuming an 
entropy of the form S{t) = A + (^gpikoslnt. Last column gives the time domain over 
which the regression is made. 



(^spiijos data points time range in hours 



0.05 0.913 ±0.002 57 40- 2,000 

0.10 0.911 ±0.008 47 100- 2,000 

0.20 0.918 ± 0.008 33 200- 2,000 

0.30 0.91 ±0.01 20 500- 2,000 



One now wants to estimate (count) the number of markers, yt{T), over a time interval 
T ending at time t. With Eq. ((Tj), this can readily be done by 



If we vary the value of t along the interval [T, NAt] , where is the total number of 
data points of the analyzed sequence. From the knowledge of yt{T), one can now readily 
calculate the probability density function p{y, T) for the moving counter. 

At this point we apply the DEA technique P^ ll2j. For each time- window of size T, 
one calculates the probability density function of yt, that is: we get p{y,T) for several 
T's. The related (Shannon) entropy, S{T), is then defined in accordance with Eq. 
In Fig. m we present the behavior of the Nord Pool entropy (using the marker defined 
in Eq. ((Zj)) for different values of threshold tq. It is observed that the entropy, S(T), 
scales in accordance with Eq. 0, and that the scaling exponent is found to be around 
^spikes = 0.91 for all thresholds (solid line in Fig. Hj). See Table [T] for further details on 
the fitting procedure. 

It is interesting to notice that a similar value (within the error bars) for the scaling 
exponent has previously been reported for the activity of the US Dollar-Deutsche 
Mark futures market (using tick-by-tick data) [25 as well as for daily index data from 
the Dow Jones Industrial Average |2n]- Hence, the Nordic power market has similar 
entropy growth to these two markets. Furthermore, all these data sets have comparable 
clustering. 

As was alluded to in the introduction to this section, the DEA technique aims at 
measuring the inhomogeneity of the distribution of the number of events (as defined by 
a marker) over a fixed period of time T, and to obtain the related scaling exponent S. 
An alternative (indirect) approach to obtaining the scaling exponent is to measure the 
tail exponent (3 of the waiting time distribution 0(r), and apply relation ^ in order to 
obtaining the scaling exponent (under the stated assumptions). 

Figure m depicts the pdf of waiting times between two consecutive spikes, 0(r), for 
different choices of the threshold tq. It is mostly consistent with a decaying power-law 
of the form (jS)). There seems to be a slight dependence of the waiting time exponent f3 



(8) 




(9) 
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Figure 4. The results of the DEA of the Nord Pool spot price data using as marker 
I ''At (01 > ^0 (with At ~ 1 hour) for different thresholds rp. The value of ^spikes is 
practically constant and, therefore, periodicities do not affect this kind of data analysis. 
For the GARCH curve, we take tq = 0.2 and we shift the results to see clearly. We 
have also added a vertical line to distinguish the intraday regime. The entropy there 
is highly affected by the consumption patterns of human activity. 



Table 2. Regressions for the tail of the waiting time pdf (/)(t) ~ l/r'^ between spikes 
larger than a certain threshold rg. The fit provides the estimated exponent /3 (see 
Eq. Qi) over the data plotted in Fig. Last column corresponds to the scaling 
exponent S predicted from the measured values of /? and assuming Eq. jSJ to be right. 





/3 


data points 


Predicted ^spikes 


0.05 


2.16 ±0.07 


39 


0.86 ±0.05 


0.10 


2.10 ±0.30 


11 


0.85 ±0.17 


0.20 


1.84 ±0.15 


11 


0.84 ±0.15 


0.30 


1.79 ±0.12 


11 


0.79 ±0.12 



on the choice of threshold, particularly for the largest waiting times r. One in Table El 
therefore presents the measured value of (3 for the thresholds used. The corresponding 
scaling exponent obtained according to Eq. (jSl), is given in the last column of Table El 
One notices that the two methods — the direct via the diffusion entropy, and the 
indirect determination from the spike waiting time distribution — are consistent within 
the error-bars. 
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X (in hours) 

Figure 5. The waiting time pdf, 0(r), vs. the time interval between spikes, t, defined 
according to |rAt(OI > (with At ^ 1 hour) for different values of the threshold ro as 
indictaed in the legend. The data collapse for the different curves towards a decaying 
power law, 0(t) ^ I/t'^, is reasonably good. The individually measured waiting time 
exponents P are given in Table [3 The solid line corresponds to a power-law fit with 
/3 = 2.16 ± 0.07 for the threshold tq = 0.05. 

5. Mean reversion and consumption patterns 

There is in tlie literature wide consensus on the fact that electricity prices are 
stationary in the sense of suffering a reverting force driving price towards a normal 
level [Ill21ElinilHll2Dlinil22l2ni- This level is most often time dependent due to marked 
consumption pattern in electricity prices caused by human social activity, and climatic 
factors like temperatures (see e.,g. Fig. We will now for consistency address the 
mean reverting character of the Nordic spot electricity prices, but now using the DEA 
technique. Previously, using the AWC techniques, the so called Hurst exponent, H, 
characterizing this property, has been measured to be about 0.4 [S], a similar value also 
found for the Californian electricity market [T]. 

For preparing for the use of the DEA for the measurement of the Hurst exponent, 
we define simply the "marker" being equal to the hourly returns, that is 



m = HPit)/P{t ~ At)) = r^t{t). 



(10) 
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The corresponding moving counting then, in analogy with Eqs. (jHl) and Q, reads 

yt{T)= ^ r^t{t) = \nP{t)-\nP{t-T) = rT{t), (11) 

i=t-T 

where one in the last transition has used the definition of the return ®. 

Figure El presents the (Shannon) entropy of the moving counter ^(t) = rAt(t). 
It is observed that the stationary state is reached after about 5, 000 hours which 
corresponds to almost seven months. The signal will diffuse anomalously until it feels 
that the motion of r/^t{t) is constrained by the reverting force. After this threshold, the 
entropy S(T) remains constant and the asymptotic scaling will he S = (cf. Eq. Q) 
and the signal loses its memory. 

According to Fig. IHl there also exists a transient state for shorter times than 
T = 5, 000 hours. In this regime, the DEA detects the presence of daily and weekly 
periodicities but it seems to be insensitive to seasonalities. The analysis, excluding the 
periodicities, also allows a fit like the one given by Eq. (jHl), that is 

S{it = A + ^return In T, (12) 

where in this case ^return equals the Hurst exponent H use to the hourly returns being 
used as markers. From Fig. IHl it is observed that the DEA calculated entropy can be 
well fitted in time by the Hurst exponent 

HbEA = ^return = 0.411 ± 0.002, (13) 

for the range between 100 and 6,000 hours. This result fully coincides with the Hurst 
exponent measured previously with the AWC method IHj. If the Hurst exponent 
H = (^return was equal to 1/2 we would talk about a (ordinary) diffusion process 
with a standard deviation growing like a/T. Therefore, finding a value smaller than 
1/2 for the Hurst exponent H = ^return means that the process is antipersistent or 
anticorrelated fI^ . 

Going further into the issue of periodicities, we now again explore the distribution 
of spikes distance r for different thresholds. We display in Fig. U\ the complementary 
cumulative probability of this random variable. Compared to probability density plotted 
at Fig. Elwe note that the complementary cumulative distribution function of the spikes 
distance is greatly distorted by the presence of periodicities in the underlying data which 
breaks a possible power law collapse for several decades. This is mainly due to the gap 
in eight hours and daily periodicities while weekly effects does not affect very much the 
cumulative curve. These distortions are consistent with the periodicities detected by 
the DEA technique. It is again difficult to see the seasonal effects but clearly detect the 
rest of periodicities. 

6. A GARCH type model 

We have seen that the spike statistics for the electricity prices is very similar to that of 
the Dow Jones index and also to the peak statistics of the US DoUar-Deutsche Mark 
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Figure 6. The result of DEA using £^{t) — rAt{t) as marker. The vertical lines 
signal the time of the main periodicities. We perform the analysis with the empirical 
data and a simulation of a GARCH model (|14(l . The two time series have the same 
growth with time although for short time lags entropy differs due to the presence of 
periodicities in the original electricity spot prices. 

futures market [211 ESI- They both have the same data clustering. The analogy leads us 
to search for good candidates for price models for electricity in a similar way to as was 
done successfully in Refs. j24| I25j. Volatility clustering and mean reversion are among 
the essential characteristics of the volatility and of the market activity. Henceforth, our 
survey focuses on existing volatility models. 

A possible new attempt is to propose a GARCH model [TH] for the spot price. We 
suggest the following model for the system spot price P(t): 

P{t) = k + xP{t-At)r]^ + uP{t-At), (14) 

where rj denotes uncorrelated gaussian random noise of zero mean and unit variance. 
Using this GARCH theoretical model, with parameters x = 4 x 10~^, u = 0.9994 and 
k = 0.00184 NOK/MWh, one can simulate long time series that results in an S{T) 
dependence that resembles that of the real Nord Pool market (Fig. 0}. To obtain the 
results of Fig. H] the total length of the simulated and real data were the same, and that 
the markers |rAt(t)| > tq were used. This model can generate Spikes — 0.9 thus giving 
similar entropy profile to the electricity spot price data. 

In addition, we also performed the DEA for the logarithmic return time series, 
generated by Eq. p4|l using the parameters given above, and markers defined in Eq. (fTIUl . 
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Figure 7. The complementary cumulative distribution of inter-spike waiting times, 
$(r) = 1 — Xxj dT'4>{T' ). One observes that the decay of these distributions are different 
in several regions. This can be attributed to the presence of periodicities in the system 
spot prices. The slope changes are important in 8 hours and day (24 hours) time lags 
but are not unobservable in a week time lag. 



Figure. El shows that this GARCH model produces results for the diffusion entropy that 
are consistent with real data. Specifically, the transient period resulting from the model 
time series has the same exponent and the stationary state is approached simultaneously 
with the real electricity return time series. Moreover, if Eq. ()14j) is rewritten in terms 
of the price difference AP(t) = P{t) - P{t - At) 

Pit - At) ^ 



APit) 



'1 



V) 



1 - V 



+ xP{t - At)r]' 



(15) 



we can observe that the reverting force has a characteristic time scale f = 1/(1 — z/). It 
is directly related to u and brings us the value f = 1, 667 hours (more than 2 months) 
which is of the same order of the time required for real data to arrive at the stationary 
state. The ratio k/{l — v) represents the normal level that prices are driven toward. 
In case, one would like to include the consumption patterns in the model one should 
replace by a periodic signal. Finally, the x value sets the magnitude of the fluctuations. 
Having a big x means having wild fluctuations of price. 
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7. Conclusions 

By applying the DEA technique we have shown that the entropy of the Nordic electricity 
spot prices grows with the size of the time window in a similar way to the volatility and 
the market activity as different as the Dow- Jones and the US Dollar-Deutsche mark 
futures. The DEA scaling parameter has the similar value ^spikes ~ 0.9 and they have a 
comparable clustering. 

We have also shown that the distribution of waiting times of the electricity price 
scales as 0(r) ~ given by Eq. and that the corresponding power-law exponent 
f3 is related to the DEA parameter ^spikes = 0.91 ± 0.01 as anticipated by the DEA - 
see Eq. (0) and Fig. El The relationship derived in Ref. |12j is obtained under the 
hypothesis that spikes are uncorrelated. This hypothesis might hold true or not but in 
any case we have found that the relationship between /? and ^spikes is still applicable in 
the Nord Pool case. 

We have also obtained a Hurst exponent (^return = 0.412 ± 0.002) with the DEA 
technique which is coincident with the one derived by previous studies and using different 
techniques. The resulting quantity implies an antipersistent behaviour of electricity 
prices. Moreover, we have been able to detect, for hourly prices, a characteristic time 
scale around 5, 000 hours in which the system reaches the stationary state. 

Due to the several similarities in statistical properties with financial markets, 
a GARCH model was proposed and investigated for the spot electricity prices. In 
particular, it was shown that the diffusion entropy of such a theoretical GARCH model 
followed surprisingly closely the entropy that was obtained for the Nord Pool system 
spot prices. 

We finally mention that although the DEA seems to detect daily and/or weekly 
periodicities in the analyzed data, a more thorough investigation is required to settle 
this issue. Efforts are currently undertaken to address this topic. 
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